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On the alleged duality of the Galactic halo 
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ABSTRACT 



We examine the kinematics of the Galactic halo based on SDSS/SEGUE data by Carollo 
et al. (2007, 2010). We find that their claims of a counter-rotating halo are the result of sub- 
stantial biases in distance estimates (of order 50%): the claimed retrograde component, which 
makes up only a tiny fraction of the entire sample, prone to contaminations, is identified as 
the tail of distance overestimates. The strong overestimates also result in a lift in the vertical 
velocity component, which explains the large altitudes those objects were claimed to reach. 
Errors are worst for the lowest metallicity stars, which explains the metal-poor nature of the 
artificial component. We also argue that measurement errors were not properly accounted for 
and that the use of Gaussian fitting on intrinsically non-Gaussian Galactic components in- 
vokes the identification of components that are distorted or even artificial. Our evaluation of 
the data leads to a revision of the estimated velocity ellipsoids and does not yield any reliable 
evidence for a counterrotating halo component. If a distinct counterrotating halo component 
exists it must be far weaker than claimed by Carollo et al. Finally we note that their revised 
analysis presented in Beers et al. (201 1) does not alleviate our main concerns. 

Key words: galaxies: haloes - stars: distances - Galaxy: solar neighbourhood - Galaxy: halo 
- Galaxy: kinematics and dynamics - Galaxy: structure 



1 INTRODUCTION 

Galactic haloes are an excellent testbed for cosmology and galac- 
tic dynamics. Their exploration can constrain the early assembly 
of galaxies as well as the dynamics of accretion of smaller galax- 
ies. Our Milky Way offers an ideal case for those investigations, 
as we can directly obtain the detailed parameters like kinematics, 
elemental abundances and physical properties of single stars sur- 
rounding us. New material is still being accreted into the Galac- 
tic halo, as the numerous streams and newly discovered dwarf 
galaxies confirm (see e.g. Ilbata et alj [l995; Kl ement et "ai] |2006l : 
iBelokurov et al.l2 007). As metallicity gradients go in lockstep with 
star formation, young accreted objects may be more metal-poor 
than old stars from the inner Galaxy and thus metallicity can not 
be simply used as a cosmic clock. This would also make it plau- 
sible that a later accreted halo component could indeed be on av- 
erage younger and more metal poor than the older parts. The for- 
mation of at least parts of the halo by accretion (combined with 
later adiabatic contraction) could give rise to differences between 
early more turbulent accretion/collapse processes and later accre- 
tion, which might leave an i mprint in difference s between the in- 
ner and outer halo (see e.g. ICooper et alJ feoiO). Another possi- 
ble source of discrepancies between inner and outer halo is dy- 
namical friction, which could be more efficient for prograde than 
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for retrograde infall dOuinn& Goodman! 1 19861 : iBvrd et aljfl986h . 
This could give rise to a different rotational signature for accreted 
material in the outer Galactic halo compared to the inner regions 
dMurante et alj2O10h . 

Historically (and as well today), halo stars have been 
extremely difficult to identify, particularly in local samples, 
e.g. demon strated by the historic arg ument between Oort and 
Stromberg dOortll 19261 : [strombergl fl927D . Like it is practically im- 
possible to get a cle an selection into thin and t hick Galactic disc 
based on kinematics dSchonrich & Binnevf2 009'), we face the anal- 
ogous problem between thick disc and the prograde stars of the 
halo. Wrong assumptions about the kinematics of the Galactic disc 
will thus affect results on the halo compone nt. Soon after the ex- 
istence of the Galactic halo was established dSchwarzschiidlll952l; 
lEggen et al.lll962h . the central question was raised if the stars of 
the inner and the outer halo had the same properties or if gradients 
or even breaks in metallicities or kinematics existed with galacto- 
centric radius. Two main strategies to identify and examine halo 
stars have been used in the past: either stars in the solar neighbour- 
hood are studied, classified according to their kinematics and then 
conclusions about the structure further aw ay are drawn by extrapo- 
lation (e.g. lSommer-Larsen & Zhenlll990h . or the surveys concen- 
trate on bright objects in the outer halo regions, s uch as RR Lyrae 
variables or globular clusters (e.g. ISandagej[T970l) . The second al- 
ternative allows to directly map the spatial struct ure by those stan- 
dard candles with good distance information (e.g. Saha 198f|). This 
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strategy implies selection biases: for example the position and pres- 
ence of RR Lyrae stars on the horizontal branch are correlated with 
metallicity and age, while it is not known if the formation of glob- 
ular clusters is representative also for all halo field stars. 

Claims of differences between the inner and outer Galactic 
halo are almos t as o ld as the discovery of the halo itself. After 
van den Berghldl967t) discussed differences in metallicity and the 
second parameter between the haloes o f the Milky Way and th ose of 
its neighbouring galaxies (M31, M33). Searle & Z inn ( 1978) found 
that Galactic clusters in the outer regions showed a larger scatter in 
the ratio of blue to red horizontal branch stars than inner halo glob- 
ular clusters, which they interpreted as a signature of an age spread. 
IPreston et alj dl99ll) found a similar difference in field BHB starsQ 
Differences in kinematics have al so been sug gested between inner 
and outer halo globular clusters dZinnlll993h . although precision 
and reliability of estimates in this respect are limited by the small 
number of available globular clusters. Various claims of an asym- 
metry in the halo azimuthal velocity distrib ution with an extende d 
tail to retrograde orbits have been made (e.g. lNorris & Rvanll 989). 
lMaiewskiUl992h even found the entire halo to be on average coun- 
terrot ating, c ould, however, not find any significant velocity gra- 
dient. iRvanl d 19921) pointed out that measurements of kinematics 
based on proper motions were particularly vulnerable to distance 
errors and showed that overestimated distances for halo stars can 

lead to false identifications of counter-rotating stars. 

In this paper we will re visit the recen t claim bv lCarollo et al.l 

d2007l) (hereafter C07) and ICarollo et alj d20ld) (hereafter CIO) 
that the Galactic halo consists of two components: a more metal- 
poor counterrotating component with larger scaleheights, distinct 
from a slightly prograde component and starts dominating the halo 
at high altitudes in their analysis In particular we will carefully 
re-examine their distance estimation procedure; we will focus on 
CIO as this paper deviates from its precursor mostly by the larger 
sample size. To avoid relying on any of the uncertain available dis- 
tance calibrations, we apply in parallel both the CIO distances and 
two native SDSS main sequence distance calibrations, checking re- 
sults additionally with an isochrone method. In Section [2] we out- 
line those methods, discuss the SDSS/SEGUE data used for this 
purpose and describe how the sample cleaning was performed. 

Thereafter (Section[3} we discuss the underlying assumptions 
and the reliability of gravity estimates used to sort stars into differ- 
ent sequences as well as the actual assumptions for absolute mag- 
nitudes by C07 and CIO. We will show that their claim to have dis- 
tances precise to 10 - 20% is unsupported and that the CIO sample 
contains a class of stars with significant distance overestimates by 

1 While it is clear that age is one parameter which will cause an older glob- 
ular cluster to be bluer than a younger one at the same metallicity, whether 
this is the dominant cause of differenc es in horizon tal b ranch morphology is 
still debated (see, for example, iDotter et afl j2010l) andfVandenBerg (2000) 
for opposing views.) 

2 This dominance of a retrograde component in the outer halo has recently 
been contested by Deason et al. (2010) although they find a retrograde mo- 
tion for metal-poor ([Fe/H] < -2) and prograde ([Fe/H] > -2) for metal- 
rich stars using a BHB sample in the outer Galactic halo. In our view, this 
issue needs to be further investigated, as they assume constant g band mag- 
nitudes for the horizontal branch in a region affected by the blue tail, which 
spans of order 2 mag in g band luminosity. If a considerable fraction of the 
halo giants is in the blue tail, their colour and temperature cuts remove a 
large part of this tail, but still leave BHB members spanning ~ 0.7 magni- 
tudes, as we tested it on SDSS photometry of metal-poor globular clusters 
known to have such a strong blue tail (M3, M13, M15 and A/92) and the 
BASTI isochrones. 



being sorted into unphysical positions in the HR diagram. In Sec- 
tion|4]we present statistical proofs of distance biases in the sample 
and in Section [5] we discuss the implications of different distance 
schemes on kinematics and the inner-outer halo dissection. 



2 THE CALIBRATION SAMPLE OF SDSS AND THE 
CAROLLO DATASET 

All the data us ed in this paper come from the Sloan Digital Sky 
Survey (SDSS. lYork eTaifc OQO), and consist of spectroscopic ob- 
servations of stars from both SDSS-I and II and from the SDSS- 
II/SEGUE survey dYannv et all 12009). Stellar parameters for the 
stars w ere estimated usi n g the SEGUE Stellar Param eters Pipeline 
(SSPP. lLee et alj2008al lb[ |Allende Prieto et alj|2008l) . 

For this study we use the calibration star sample from SDSS 
public data release 7 (DR7)0As the colour transformations and dis- 
tances used in C07 and C10 are not part of the public data releases, 
we draw this information from the sample used by C10, which is 
a cleaned version of DR7 and which was kindly provided on our 
request. 

The calibration star sample comprises two datasets, namely 
the photometric calibration star sample and the reddening calibra- 
tion star sample. Querying the DR7 catalogue for these stars pro- 
duces a total of 42841 entries, but many refer to identical objects, so 
that the actual number of unique objects in the database is 33023. 
The interested reader is referred to the Appendix where we describe 
how the cleaning of the sample from questionable objects has been 
performed that cuts down the sample to 28844 stars. 

Throughout the paper we make use of two different classes 
of distance determinations (see Appendix for details): On the one 
hand the distances used b y C07 and C10 and on the other hand two 
derivations adopted from llvezic et al.l (2008). When using the dis- 
tances by C07 and C10 we examine the effects of their different 
sequences by both using their entire sample ("Carollo all") and us- 
ing exclusively their dwarf stars ("Carollo dwarfs"). For the native 
SDSS calibrations we restrict ourselves to the dwarf stars, impos- 
ing in general a gravity limit of log(g) > 4. 1 to reduce the impact 
of giant and subgiant contamination). We there use two different 
schemes: An adopted main sequence derivation, where we increase 
distances defined via their Eq. (A2) and (A3) by accounting for al- 
pha enhancement and additionally decreasing the adopted absolute 
magnitudes of all stars by 0. 1 mag, hereafter termed the adopted 
main seque nce calibration, sho rt "IvzMS"). Second the calibration 
favoured bv llvezic et al. I d2008l) using Eq. (A7) from their appendix, 
which is steepened towards the isochrones in order to account for 
age dependent effects, hereafter called the age dependent calibra- 
tion (short: "IvzA7"). We would like to point out that we believe 
neither distance calibration to deliver the full truth, yet they are cur- 
rently the most commonly used schemes and can give hints on the 
intrinsic uncertainties of all methods. To have an additional test we 
cross-checked and confirmed our findings with distances derived 
directly fro m using 12.5 Gyr BASTI isochrones dPietrinferni et al.l 
l2004 l2006); some details are provided in the Appendix. 

C10 applied two different geometric cuts. The first was a dis- 
tance cut at an estimated distance of 4kpc, limiting the errors on 
velocities caused by proper motion errors. The second was a cut 
that removed stars differing from the Sun by more than 1.5 kpc 
in galactocentric radius (i.e. outside 7.0 kpc < R G < 10 kpc with 

3 http://casjobs.sdss.org/dr7/en/ 



On the alleged halo duality 3 




150 200 250 300 350 
Galactic I 

Figure 1. The distribution of all reliable stars in the calibration sample, 
fulfilling the temperature range limit and with acceptable kinematic infor- 
mation. To demonstrate the weak line-of-sight motion support of azimuthal 
velocities we plot the 0.95 (innermost lenses, green), 0.8, 0.5 and 0.3 (out- 
ermost lenses, light blue), contours of the angle term r]y in Eq.jTJ that quan- 
tifies the relative support by the direct line-of-sight velocity measurements. 
This demonstrates the vulnerability of V velocities in this sample to distance 
errors as they cannot be found directly from line of sight velocities. 



their value of R a& = 8.5 kpc for the Sun). The latter cut was rea- 
soned in C10 by the use of their ap plied orbit model, which was 
adopt ed from Chib a & Beers (2000) using a Stiickel-type potential 
from[s ommer-Larsen & Zhenl a 1990I) . As we do not make use of 
this orbit calculation, in this work we do not apply the second cut, 
which removes of order one third of the available stars. 

Due to the magnitude ranges in the SEGUE survey the dis- 
tance cuts imposed by C10 and C07 remove almost all giants from 
the sample, a minor fraction of subgiants and very few dwarf stars. 
The effects of both selections depend strongly on the adopted dis- 
tances and thus the absolute magnitudes of the stars. Applying the 
C10 distances to the sample, 21600 stars remain within the distance 
limit of 4 kpc, of which 14763 have surface gravities log(g) > 4.0 
and are thus classified as dwarf stars. With the adopted main se- 
quence calibration 15808 stars with log(g) > 4.0 are found in the 
sample, which gives a first hint to the more stretched distance scale 
by C10. 

Metallicities are taken from the DR7 pipeline adopted values 
dLee et all l2008al fbh- For a discussion of the different metallicity 
scales of DR7 and C10 the reader is referred to the Appendix. There 
we also describe in more detail how kinematics are derived from 
distances, radial velocities and proper motions. 



3.1 Effects of distance errors 



Selecting a star into the wrong position in the colour-magnitude 
diagram results in a faulty estimate of its absolute magnitude and 
thus an erroneous distance. As they are the most common popula- 
tion we would naively expect the largest contamination to be main 
sequence stars mistakenly addressed as turn-off stars. These will be 
assumed to be far brighter than they are, hence their distance will be 
overestimated, bringing many of them, especially halo stars, falsely 
into the retrograde tail of the velocity distribution. This effect hap- 
pens via the translation of proper motions into velocities and thus 
prevails for samples that have low support by radial velocities (cf. 
Section [8T1 of the Appendix). Fig. Q] depicts the locations of stars 
(red dots) in Galactic longitude (x-axis) and latitude (y-axis). The 
SDSS/SEGUE sample is (due to the location of the telescope in the 
northern hemisphere and the strategy to avoid the high extinction in 
the Galactic plane) largely concentrated away from the plane and 
towards the Galactic North Pole. Consequently it has almost no 
points at directions that would have high support of azimuthal ve- 
locities by the direct line-of-sight velocity measurements. The con- 
tours in the plot encircle the regions of high line-of-sight velocity 
support of V/, in the sky0 Within the ellipses, the fraction rj v (Eq.lJTJ 
in the Appendix) of the line-of-sight velocity going into the helio- 
centric azimuthal velocity V;, is larger than 0.95 (smallest lenses), 
0.8, 0.5 and 0.3 (largest lenses). They demonstrate how heavily any 
analysis of the azimuthal velocities has to rely on the transverse ve- 
locity component and thus proper motions and distance estimates. 



The effect of distance errors in this process is easily under- 
stood: Think of driving a car past a field that has a rabbit sitting on 
it. As the speed of the car is known, the fact that the rabbit rests on 
the lawn can be derived by the car driver from its apparent angu- 
lar speed - if the distance is right. If the natural size of the rabbit 
is overestimated, so will be its distance and to explain its angular 
motion one wrongly infers that it moves opposite to the car's direc- 
tion of motion. And vice versa a distance underestimate drags the 
estimated rabbit velocity towards that of the car, i.e. it is wrongly 
inferred that the rabbit moves in the same direction as the car does. 
As the Sun moves with a velocity of more than 200kms~ 1 around 
the Galactic centre, for this sample a 10% distance overestimate im- 
plies that an average halo star in the sample will be wrongly pushed 
by ~ 20 km s -1 into retrograde motion, and a larger distance error 
will entail a proportionately larger retrograde motion. 



3 ASSESSING THE DISTANCE CALIBRATIONS 

The core assumptions in C10 are those about stellar distances. They 
claim that by using the log(g) estimates from the SEGUE stellar 
parameter pipeline stars can be reliably sorted into clean sequences, 
i.e. main sequence , turn-off and subg iant/giant. The attraction of 
this idea, laid out in iBeersetal J ( 120001) . is that it seems to reduce the 
distance errors to simple uncertainties in colour and metallicity on 
well-determined sequences. However, as we will see below, things 
are more complicated. 



It should also be mentioned that azimut h al vel ocities are sub- 
ject to an error similar to the iLutz-Kelkerl d 19731) bias: Even if 
there is no net bias on estimated absolute magnitudes and thus of 
distances, a symmetric magnitude error distribution will cause an 
asymmetric distribution of estimated distances with a longer tail 
in the overestimates. As these directly translate to the transverse 
velocities one will with any magnitude based distance scheme en- 
counter asymmetric velocity errors that give rise to an extended 
counterrotating tail in the measured halo azimuthal velocity dis- 
tribution. This may a lso be related to the asymmetry found by 
iNorris & Rvar] Jl989l) . 
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Figure 2. The distribution of the "well behaved" stars in the calibration 
sample (red crosses, without distance cut to make the giants visible) in the 
dereddened colour-gravity (DR7 pipeline) plane. Horizontal lines mark the 
cuts adopted by Carollo et al. between the dwarf, turn-off and subgiant/giant 
regions. Stars that have likely membership in their counter-rotating halo, i.e. 
with Vi, < -400 km s and fulfilling d < 4kpc are marked with blue points. 

3.2 Gravities 

Fig.[2]shows the positions of all cleaned sample stars in the colour 
- log(g) plane (red points). The horizontal lines show the selection 
regions in log(g) used by C10. Stars with log(g) > 4.0 are classi- 
fied as main sequence stars. Those with log(g) < 3.5 are assumed 
to be subgiants or giants, while the intermediate objects are clas- 
sified into their "turn-off' branch. The blue dots mark the objects 
in the C10 sample that have heliocentric azimuthal velocities, V/,, 
smaller than -400kms~', i.e. they are on highly retrograde orbits 
and should be in the majority counter-rotating halo objects accord- 
ing to C10. There can be minor differences between the gravities 
directly from DR7 and their sample, as we we use the best deter- 
mined values when their are double or multiple entries for a star, 
while their log(g) values (which we do not have access to) are prob- 
ably averaged. 

The plot reveals one crucial problem with those criteria: al- 
though it may be expected that some distinct branches of stars are 
present, they are not reflected by the SEGUE measurements. While 
the upper giant branch is apparent, we cannot make out a substan- 
tial decrease of densities between the main sequence and subgiant 
regions. Measurement errors prevail, especially in the turn-off re- 
gion, and veto against a clean selection of the components. 

One possibility is that the measurement errors on log(g) are so 
large that it is not possible to use the log(g) measurements to clas- 
sify stars into these three categories. Also the true main sequence 
is inclined in gravity versus colour. In this perspective the constant 
(colour and metallicity independent) gravity cuts applied by C07 
and C10 do not appear well founded. If one aims to select a pure 
dwarf sample, the tightening of the constant cut relative to the in- 
clined main sequence, is, however, beneficial in reducing the turn- 
off-contamination. 

The accurac y of the DR7 log(g) values is discussed by 
iLee et aD {2008b), who show estimated surface gravities for open 

4 Throughout the paper we distinguish between the velocities in the solar 
rest frame and coordinate system (C/,, V/,, W;,) and velocities in the Galactic 
rest frame and Galactic cylindric coordinate system (U, V, W). For a short 
discussion of those we refer to the Appendix. 



and globular cluster stars observed by SEGUE f or the pu r pose of 
calibration of the survey. As demonstrated by ILee et al.l ([2008b ) 
(see their Fig. 15 ff.) the low spectral resolution results in a signifi- 
cant scatter in log(g) values: some turn-off stars have log(g) < 3.5, 
and significant numbers of stars with log(g) measurements in the 
region assigned by Carollo et al. to the turn-off (3.5 < log(g) < 4) 
are clearly either subgiants or main sequence stars. Thus we have to 
expect a sig nificant number of main se quence objects in the turn-off 
band of the Carol lo et al.l d2007l,l2010h samples. 

The consequent bias in distance estimation discussed above 
will give an artificially enhanced fraction of retrograde halo stars 
residing in the turn-off/subgiant regime. Although Fig. [2] should 
play this effect down with the red points being drawn from the 
main sample without distance cuts (to keep the red giant branch 
visible), while the blue circles satisfy d < 4kpc in addition to the 
velocity cut, this crowding of the counterrotating halo stars into the 
designated turn-off region is still prominent. The distance cut is re- 
sponsible for the missing enhancement in the giants. 

Carollo et al. also find that their outer counterrotating halo 
members display significantly lower metallicities than the average 
of the inner halo. Could this be related to problems with distance 
estimation as well? It would not be unreasonable to expect the ac- 
curacy of the log(g) estimates to decrease for the stars of lowest 
metallicity, increasing by this the fraction of fake turn-off stars and 
subgiants, since the stellar lines become weaker. In fact, Ma et 
al. (in prep.) find this effect. A better understanding of why the 
"counter-rotating" component of Carollo et al. has lower metal- 
licity can be achieved from Fig. [3] In its top panel we plot the 
gravities against the metallicity, again with red points for the entire 
sample and with blue circles for the strongly retrograde stars. Stars 
with lower metallicities in the sample get on average assigned lower 
gravities. This could arise from the fact that determining gravities 
gets more difficult on the lowest metallicity side. The inclination of 
the sample in the metallicity-gravity plane favours distance overes- 
timates for metal poor stars, thus lowering the average metallicity 
of the "counter-rotating" component. This shift in classification of 
stars is demonstrated in the middle panel of Fig.[3] Almost all iden- 
tified subgiants are metal-poor, while the average metallicity rises 
towards the higher surface gravity categories. The bottom panel 
shows the ratio of the number of turn-off to main sequence stars, 
p, against metallicity. One might argue that the pronounced rise 
of p towards lower metallicities be caused by the intrinsic rarity 
of metal-poor objects. As the sample is dominated by magnitude 
cuts, the different ratios could, however, just be explained by dif- 
ferent geometry of the subpopulations, which indeed puts stars with 
higher scaleheights (e.g. halo) to more remote positions than those 
with low scaleheights (e.g. disc). Since the C10 sample extends by 
definition less than 4 kpc away from the plane it is hardly possible 

to explain how p can rise by a factor of ~ 3 from \FejH\ D! a 1.5 

to [Fe/H] DR7 ~ -2.5 as the large scaleheight of the halo should 
veto against relative density variations of the halo populations by 
this amount. 

3.3 Absolute magnitudes 

Having discussed the surface gravity estimates it is time to turn to- 
wards the actual assumptions on absolute magnitudes of stars from 
which the distances are inferred. The claim by C10 to reach an ac- 
curacy of "10 - 20%" in distance estimates is predicated on the 
ability to cleanly select stars according to their spectroscopically 
determined gravity into several branches. If such a selection was 
feasible, it would indeed limit the uncertainties to those caused by 
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Figure 3. Top panel: The distribution of the cleaned "all star" sample in the 
metallicity- surface gravity plane. Horizontal lines mark the cuts adopted 
by Carollo et al. between the dwarf, "turn-off" and subgiant/giant regions. 
Stars in the CIO subsample displaying a clear membership to what might be 
a counter-rotating halo, i.e. with Vi, < -400.0 km s are marked with blue 
points. Middle: Metallicity distributions for the different gravity classes 
used by Carollo et al., showing how low metallicities dominate the "turn- 
off" and "subgiant" stars. Bottom panel: Ratio of identified "turn-off" stars 
to dwarf stars as function of [Fe/H^oRl using the gravity cuts by Carollo 

et al. Even between [Fe/H]oR7 1-5 and [FelH]oRl 2.5 there is a 

strong uptrend. 
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Figure 4. The distribution of the sample stars in the (g - i)o,M r plane 
with tiny crosses. Colours give the metallicity [Fe/H]c a r according to the 
colour bar on the right. The top panel contains all stars in the C10 sample 
that pass our quality cuts, while the middle panel exclusively contains stars 
with [FelH]Dia < -1.9. For comparison we give the BASTI (Pietrinferni et 
al. 2004, 2006) alpha enhanced isochrones at [Fe/H] = -2.14 of ages 10 and 
12.5 Gyrs (green lines) and the Ivezic et al.(2008) absolute magnitude cali- 
bration at the same iron abundance (red line). In the bottom panel the metal- 
poor stars have been enhanced in size and stars with V;, < -400 kms -1 de- 
rived with C10 distances are marked by light blue circles. Red lines show 
the adopted main sequence calibration (IvzMS) at [Fe/H] = -2.14 and 
[Fe/H] = 0. Dark blue lines show the BASTI isochrones at solar metal- 
licities for 1, 4, 10 and 12.5 Gyrs, green lines same as in the middle panel. 
The crosses mark the typical photometric errors and uncertainties in the 
isochrones. 
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metallicity, reddening and photometric errors. The reader should 
bear in mind that the low relative density of the outer halo compo- 
nent found by Carollo et al. implies that even a contamination on 
the 1% level (i.e. 200 out of ~ 20000 stars) can alter the results. 

The sorting into different branches via the formalism of Beers 
et al. (2000) can be seen in the top panel of Fig. [4] The plot shows 
all stars in the C10 sample in the (g - i)o, M r plane that do not have 
warning flags. The absolute magnitudes were derived directly from 
the C10 distances using the distance modulus via 



M r = 7- -51og 10 ( 



0.01 kpc 



(1) 



where M r is the derived absolute magnitude, r is the reddening 
corrected apparent magnitude from the DR7 database and dc ar the 
distan ce given by Carollo et al. and derived from the iBeers et al.l 
J2000l> sequences. 

The three branches, into which the sample stars are selected, 
can be identified in the figure. Lower metallicity stars are brighter 
and bluer, but the latter shift dominates. Thus the main sequence 
gets fainter at the same colour. As the metallicity spread of es- 
pecially the main sequence in the top panel of Fig. [4] partially 
obscures the underlying sequences, we restricted the sample to 
[Fe/H] D/ n < -1.9 in the middle panel. For comparison we plot 
the adopted main sequence calib ration at [Fe/H] ~ -2.14 and the 
alpha enhanced B asti isochrones dPietrinferni et a l. 2004, 20Q(J) for 
SDSS colours (cf . iMarconi et aLlbood) at TFe/Hl = -2.14 at ages 
10 and 12.5 Gyrs. The main sequence and the giant branch are 
apparent and between them lies a strong sequence sloping down 
from (0.2, 3.8) to about (0.8,4.8). The stars in this band are termed 
turn-off stars by Carollo et al. as they comprise everything that has 
gravities 3.5 < log(g) < 4. As we can see from the bottom panel, 
this artificial turn-off sequence comprises the majority of heavily 
counter-rotating stars (green dots). The reader might ask why we 
do not see a large number of stars as false identifications in the 
brightest subgiant branch. As already discussed above the abso- 
lute magnitude difference to those stars is, however, so large that 
by the distance limits most stars sorted into the subgiant and giant 
branches will be dropped from the sample, as they are deemed to 
be more than 4 kpc away. This conclusion is verified by the sub- 
giant branch getting much more populated when the distance cut is 
removed (cf. Fig 2). 

In stellar evolution there exists no unique turn-off branch. In- 
stead there is just a region where stars leave the main sequence and 
move up to the subgiant branch, and which depends on metallic- 
ity and age of the population. Thus the artificial sequence we face 
here can only be thought of as a compromise for stars in the broad 
region between main sequence and subgiant branch that intends to 
describe the average luminosity of these objects. In this perspective, 
the claim by Carollo et al. to achieve 10 - 20% accuracy in dis- 
tances is ruled out in this transition region as it intrinsically spans 
more than one magnitude (at fixed metallicity and colour), i.e a dis- 
tance uncertainty of more than 50%. Because this "turn-off' branch 
is constructed to be a compromise between subgiants and main se- 
quence, the gap between the subgiant branch (moving up towards 
lower metallicity) and the main sequence (moving down, as the 
colour effect dominates) widens with decreasing metallicity. This 
aggravates the effects of misassignments on estimated distances for 
lower metallicity objects. 

The worst problem with the "turn-off" stars appears in com- 
parison with stellar models. Comparison of the adopted positions 
of metal-poor stars in the colour-magnitude diagram with the 
isochrones at 10 and 12.5 Gyrs reveals that most objects in the arti- 



ficial "turn-off" branch reside on the red side of the low metallicity 
turn-off region. Some of the selected "tum-off" stars are even on 
the red side of what could be achieved at solar metallicity. In fact 
the overwhelming majority of strongly retrograde stars is actually 
claimed to be in a region where according to our knowledge of as- 
trophysics no star of reasonable ages can reside. 

3.4 Colour transformation and main sequence comparison 

It should be mentioned that the distances used in C07 and C10 
were derived using the B - V co l our ca libration and thus the 
colour transformation of iLee et alj d2008al) had to be performed 
to apply them to SDSS colours. F ortunately there are now good 
isochr ones (e.g . BAS TI isochrones, IPietrinferni et alj |2004). fidu- 
cials dAn et al.l l2008h and photometry in SDSS colours readily 
available, making such a colour transform ation to tran s late SDSS 
colours to the former B - V calibration bv lBeers et alj d2000h un- 
necessary. The colour transformation may explain some of the 
scatter and systematic shifts presented in Fig. [5] that shows the 
comparisons of the absolute magnitudes from C10 to the native 
SDSS calibrations. The latter is depicted by red lines in Fig. [4] at 
[Fe/H] = -2.14 and [Fe/H] = 0. Both isochrones and the two main 
sequence approximations are fainter than the adopted C10 abso- 
lute magnitudes, especially at lower metallicities. This can be seen 
from the middle panel of Fig. [4] and from Fig. [5] which depicts the 
distribution of differences in absolute magnitudes of C10 towards 
those derived via the adopte d main sequence ca libration (top pan- 
els) and those derived via the llvezic et alj {2008) age-dependent re- 
lation (bottom panels) for all stars (red), metal poor objects with 
[Fe/H] < -1.9 (blue) and for the "counter-rotating" stars with 
V h < -400kms~' (green crosses and errorbars). Errorbars depict 
the Poisson noise. On the right side of Fig. [5] we show the same 
quantities, but ex clusively for redder stars with (g - ;)o > 0.4 where 
the sequences of IBeers et all (2000) have a larger separation, and 
which also excludes the expected turn-off region (cf. the isochrones 
in Fig. [4} for the very metal-poor stars. Apart from the observa- 
tional scatter, the peaks arising from the main sequence A M ~ 0.3 
and from the "turn-off branch" around A M ~ 1.0 are washed out for 
the whole population, as the sequences shift with metallicity and 
the main sequence dominates. However, the generally higher intrin- 
sic brightness under the Carollo et al. assumptions is clearly seen 
in both distan ce descriptions. On the blue end the age-dependent 
formulation bv llvezic"e t al. (2008) is a bit brighter than the adopted 
main sequence calibration and makes the peak of the "turn-off se- 
quence" in the metal poor and counterrotating subsamples of C10 
merge with the main sequence at blue colours, yet the large offset 
remains. 

For stars that are claimed by C10 to have heliocentric veloc- 
ities V h < -400kms~' and thus make up the bulk of what they 
identify as a counter-rotating halo the two peaks of main sequence 
and turn-off stars are separated, as the metallicities of these stars 
are more narrowly distributed . For comparison we show the result 
of selecting only stars that have [Fe/H] DR j < -1.9 (blue squares 
in Fig. [5}. These distributions look very similar apart from the 
counter-rotating stars being more inclined to the bright side. Also 
regard the increased importance of the turn-off band for those pop- 
ulations. 

In summary there are two main drivers of overestimated dis- 
tances in Carollo et al. compared to the native SDSS distance cali- 
brations: The main sequence of C10 for metal-poor stars is brighter 
than in the other calibrations and the number of stars in the turn- 
off and subgiant/giant branches is increased by a large factor for 
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Figure 5. The distribution of absolute magnitude differences between the Carollo et al. calibration and the other metallicity-dependent colour-magnitude 
calibrations. In the top row we show the results for the adopted main sequence calibration ("IvzMS"), while in the bottom row we compare to the age- 
dependent calibration by Ivezic et al. (2008) ("IvzA7"). On the left we show the entire colour range, while we plot only redder objects with (g - i)o > 0.4 on 
the right. In the entire sample (red bars) the main sequence by Carollo et al. is brighter by ~ 0.15 mag than in the native SDSS calibrations. The offset gets 
larger for the counter-rotating (green crosses) and metal-poor (blue squares) stars; note that the Carollo et al. turn-off sequence is very prominent with an offset 
of about 1 magnitude. 



metal-poor stars; most of the designated "turn-off" stars by C10 re- 
side in positions in the colour-magnitude plane that would require 
unreasonable ages. We thus conjecture that the "counter-rotating" 
halo as presented by C10 is due to distance uncertainties and the 
selection of stars into unphysical stellar branches an artifact. 



4 SIGNATURES OF MISSELECTION 
4.1 Geometric distribution anomaly 

Even if the biases in distance assignments were not obvious there 
are still ways to detect them. As erroneous distance estimates 
mostly act on the part of motion perpendicular to the line of sight, 
populations with distance errors will be preferentially found in re- 
gions of the sky perpendicular to the velocity component in ques- 
tion, or in other words, if stars end up in one component by miss- 
election this should should show up as a bias in sample geometry. 
Such a bias will not be aligned directly with the Galactic coordi- 
nates, but with the part of the biased velocity component covered 
by proper motion. This is obvious as the largest errors should hap- 
pen where the uncertainty in the motion is largest, i.e. where the 
proper motion and distance estimates have the largest impact. For 




?pm, v 

Figure 6. The cumulative distribution of stars on the part of Vi, that is cov- 
ered by proper motions and not radial velocity measurements. The stars 
with Vi, < -400 km s are drawn with a purple line while the overall sam- 
ple distribution is drawn with a red line. Beneath is the 0.01 significance 
level of a Kolmogorov-Smirnov test. 
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these statistics we make use of the squared angle terms connecting 
the proper motion to the azimuthal velocity (from now called the 
"proper motion partition"). We take from eq.© in the Appendix 
the angle terms that connect V A to the proper motions / and b, square 
them and add them together: 

(„„ hV = sin 2 (/)sin 2 (6) + cos 2 (/). (2) 

Fig. [6] shows the cumulative distribution of stars over the proper 
motion partition for V h (red line). The CIO sample is (as the en- 
tire SEGUE survey) concentrated towards positions where the az- 
imuthal velocity is mostly covered by proper motions as it is mostly 
oriented towards high Galactic latitudes. Yet the counter- rotating 
subsample (purple line) is even more concentrated towards the high 
proper motion contributions. The average value for ( pm ,v rises from 
0.909 for all stars to 0.933 for the subsample with V h < -400.0. In 
other words, the fraction covered by robust radial velocities drops 
from 0.091 to 0.067. In a Kolmogorov-Smirnov test the probability 
for equality of the two distributions is well below the 1 % level and 
thus equality is strongly rejected. Again the high values for £, pm y 
show how vulnerable the sample azimuthal velocities are to any 
distance errors. 

This is of course not a proof of the bias, but a strong indica- 
tion. Against this argumentation one could raise the objection that 
stars at low metallicities are located more polewards, as the halo to 
disc ratio in the sample rises (higher altitudes are reached) and so 
metal-poor / halo stars show a different spatial distribution. Indeed 
for lower metallicities the sample distribution shifts polewards re- 
ducing the difference, which remains present in all cases, but can 
become insignificant due to the shrinking sample numbers. 

4.2 Linear error analysis and velocity crossterms 

A robust approach to prove and quantify the distance errors 
is presented in Fig. [7] which shows the W>, velocities against 
sin(Z) s'm(b) cos(i) for all stars that should be part of the counter- 
rotating halo, i.e. have V h < -400kms~'. The strong uptrend can- 
not be any stream, which would show up as a narrower band. It 
appears that the entire halo experiences a lift in W velocity, which 
has its origin in the distance errors of the C10 analysis. By wrongly 
assessing distances the relative motion of stars in the azimuthal di- 
rection mixes over the proper motions into the vertical component 
depending on this angle term. To illustrate this effect think of a star 
that is seen at b = 45", / = 90° and has U h = 0, V h = -500 km s" 1 , 
Wh = 0. Due to the geometry of the setup the star will have a sig- 
nificant proper motion in b. With a wrong distance estimate the 
relative impact of line-of-sight velocity and proper motion changes 
and the star is assigned a non-zero W velocity. A more detailed 
description is found in Schonrich et al. (in prep.). The fitting line 
has a slope of 119.5 + 16.5 kms~', i.e. the trend is significant at a 
level of more than la. We can thus state that a significant distance 
error that is expected from the discussions above is here quantita- 
tively demonstrated. Since the alleged counterrotating halo stars in 
this subsample (with a cut at Vh < -400kms~') have on average 

V], 470kms~ 1 , we arrive at a distance bias of roughly 40% for 

the counter-rotating stars, which translates into a magnitude error 
of around 0.7 mag. It should be mentioned that errors might even be 
higher: Selecting more strongly retrograde stars the trend estimate 
gets even higher. 

Restricting the sample of stars with V;, < -400kms~' further 
to dwarfs (log(g) > 4.0) and thus removing the largest identified 
source of distance errors, the numbers of stars in the sample drop 
by the selection dramatically from 735 stars to 299. The removal of 
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Figure 7. The counter-rotating stars display an intense trend in average W/, 
velocity over the angle combination that describes the connection of the 
azimuthal and vertical velocity components via distance errors. 



the low gravity stars also gives fewer outliers in the velocity distri- 
bution. Only one object has \W\ > 400 km s -1 . The lower panel of 
Fig. [7] shows the remaining stars together with an equivalent linear 
fit. Again the trend is highly significant, though now more moderate 
with 83.8±22.5kms-'. The <t w in the sample is already down from 
129 km s -1 to 110 km s~ ] . The more moderate trend corresponds to 
an error of about 25%. It also fits well into the picture that the dis- 
tance overestimate for the full sample is larger than for the dwarfs, 
when we remove the spurious "turn-off' branch. This bias has the 
effect that the distance overestimates do not only increase the W 
velocity directly, but also give rise to an additional contribution to 
the W velocity dispersion by turning a part of the large heliocentric 
azimuthal motion of halo stars into a fake vertical term. We will call 
this behaviour "velocity crossovers". This effect is even more im- 
portant in the original analysis by C07 and C10, as they restrict their 
sample in galactocentric radius, thus increasing the weight around 
/ = 90° and I = 270°, where sin(Z) is largest. 



5 VELOCITY DISTRIBUTION 

How do the described biases and different distance derivations af- 
fect the velocity distribution and what might remain of the alleged 
counter-rotating component halo when we reduce the distance bi- 
ases that spuriously inflate it? Fig. [8] shows the velocity distribu- 
tions under different cuts and distance approximation both on a 
linear scale and on a logarithmic scale that reveals the wings of 
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Figure 8. Velocity distributions for stars in the calibration sample. To make 
them comparable at the different sample sizes, all distributions were nor- 
malized to unity (i.e. divided by the total number of stars in each subsam- 
ple). Error bars show again Poisson errors. 



the distributions. Error bars give the Poisson noise, but neglect 
any other sources of error. The velocity distributions of the halo 
and the Galactic disc can be made out in both plots. The long 
tail of the counter-rotating halo is clearly visible for the C10 sam- 
ple (red circles). As already noted above, a restriction to dwarf 
stars (green crosses) diminishes the counter-rotating tail (around 
V < -170kms~', which corresponds to V h < -400kms~'). The 
blue filled squares show what we obtain with the adopted main 
sequence calibration, and the purp le empty circles dep ict the age- 
dependent distance calibration by llvezic et"al] | |2008|) . With both 
native SDSS distance calibrations the sample displays a clean 
downtrend to the low velocity side corresponding to a V velocity 
dispersion of around 70 - 80 km s -1 . 

One might argue that the two native SDSS calibrations lead to 
some contamination by giants that are treated as dwarfs and thus 
provoke the opposite effect of underestimating distances. To re- 
duce the contamination by evolved stars we use a gravity cut of 
log(g) > 4.1 in the following discussions, which mostly helps to 
reduce the density saddle between the disc and halo component, 
where halo stars with severe distance underestimates tend to assem- 
ble. However, the left tail of the halo distribution does not change 
significantly on a tightening or loosening of the gravity cut. In the 
following discussions we will always use the slightly tightened con- 
dition log(g) > 4.1, but we checked that all of our conclusions are 
valid regardless of the specific choice of the gravity cut. The ab- 
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Figure 9. Vertical velocity dispersions against metallicity for different halo 
subsamples and distance calibrations. In all samples we cut away all pro- 
grade stars to eliminate most of the disc contamination. While the Carollo 
et al. full sample (red line) shows a vigorous uptrend towards the lowest 
metallicities, this trend almost vanishes in their dwarf star subsample (green 
line), while using the two native SDSS calibrations no significant trend is 
detectable. 



sence of a significant excessive tail for dwarf stars can in our view 
only be explained by the fact that strongly counter-rotating halo 
dwarfs are at the best an extremely rare population. It is highly 
implausible that the counter-rotating halo consists exclusively of 
subgiants and giants. 



5.1 Kinematics versus metallicities 

The reason why Carollo et al. (2007, 2010) get a transition be- 
tween their inner and outer halo from a local sample is the in- 
creasing scaleheight of their populations with lower metallicity, 
which is linked to the increasing velocity dispersion perpendicular 
to the plane. According to our above discussion of kinematic finger- 
prints of distance errors, at least the increase of vertical dispersion 
by velocity crossovers (Section [4.2| l due to distance overestimates 
should disappear when cutting away the unphysical turn-off stars 
and further when switching from C10 dwarf distances to the native 
SDSS calibrations. This is indeed observed in Fig. [9] It shows the 
values of vertical velocity dispersion against metallicity for differ- 
ent subsamples together with Poisson errors. The full C10 sample 
(red circles) harbours a prominent uptrend towar ds lower metal- 
licity, which is at odds with the earlier result from Chiba & Beerj 
(2000). This uptrend almost vanishes in their dwarf subsample 
(green crosses). Both when using the adopted main sequence cali- 
bration (blue filled squares) and when using the llvezic et a l] d2008t) 
age-dependent distance estimates (purple empty squares), we see 
no significant trend with metallicity any more; this finding is robust 
against changes in the gravity cut. To guide the eye we plot a hor- 
izontal line at 85 km s . There is a suggestion of a subtle increase 
below [Fe/H] DR1 < -1.8, yet numbers are too small to allow for a 
judgement. A source of uncertainty is the degree of possible giant 
contamination. It likely rises towards the metal-poor side, which 
may result in a reduction of measured dispersions. On the other 
hand the radial velocity support for the vertical velocity compo- 
nent is very high due to the polewards orientation of the sample. So 
the decrease effected by distance underestimates is relatively weak. 
Underestimating distances takes a further moderated effect on ver- 
tical dispersions: While a distance overestimate both increases the 
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measured vertical dispersion by direct overestimate of the proper 
motion part, and a lift via the velocity cross-overs, on distance un- 
derestimates the effects have opposite signs, thus to some part bal- 
ancing each other. 

Another way to look at the problem is to plot metallicity dis- 
tributions as a function of the kinetic energies (by stellar mass) of 
stars. In Fig.QJJwe show the separations in entire kinetic energy 
(v 2 , left column) and in vertical kinetic energy (W 2 , middle col- 
umn). In the rightmost column we show the vertical kinetic en- 
ergy, but for retrograde stars (V < -10kms~'). As the rotation 
velocity of disc stars adds to the kinetic energy the disc has its 
strongest dominance of course not in the lowest energy bin, but 
at 180 2 (kms~') 2 < v 2 < 270 2 (kms -1 ) 2 , as most disc stars are on 
quite circular orbits with an entire velocity close to the circular ro- 
tation speed of the Galaxy. The prominent shift especially in the 
left tail and also peak of the metallicity distributions in the orig- 
inal CIO distance prescription for the full sample (top row) can 
be seen in both the total and vertical energy separation. However, 
this already diminishes, when we plot the subsample of the CIO 
dwarf stars with log(g) > 4.0 (second row). No convincing outer 
halo signature can be found in either of the native SDSS calibra- 
tions (bottom row). At the highest vertical energies a single spike 
at [Fe/H] ~ -2.2 sticks out. At first glance this could be taken as 
a hint for an outer halo. However, in this interpretation it should 
be mirrored by the entire energy distribution, which is not the case. 
We are tempted to identify this feature at l east partly wi t h a pr omi- 
nent metal-poor stel lar stream described bv lHelmi et al J ( 119991) and 
Keplev et al] d2007l) . which we expect exactly at these very high 
vertical energies. This stream can as well serve as explanation for 
the subtle and insignificant increase of vertical velocity dispersions 
at [Fe/H] DI a < -1.8 that was seen in Fig.|9] 

There appears also a slight general drift towards lower metal- 
licity at the highest energies, which can be mainly seen in the verti- 
cal term. Caution should be exercised, however, due to the varying 
presence of the Galactic disc (most prominent at the lowest verti- 
cal kinetic energies and at total energies corresponding to the ro- 
tation speed), which then impacts on the normalization of the halo 
component and fools the eye because the apparent halo peak can be 
shifted by the wing of the disc distribution. A quite robust approach 
is removing all prograde stars from the sample to minimize disc 
contamination. As this biases kinematics, the overall energy distri- 
butions are altered, but the vertical energies (right hand column) 
should not be affected. In Fig. [TO] the entire CIO sample shows 
a clear signature of lower metallicities in the higher vertical en- 
ergy bins. When removing the contested turn-off stars (i.e. plotting 
dwarfs, centre row) no trend apart from the discussed spike can be 
detected regardless of the applied distance determination. 

In summary it can be stated that also in terms of metallicity 
versus energy the sample has no reliable outer halo signature, nei- 
ther when using the native SDSS calibrations nor when using the 
CIO dwarf subsample. 

5.2 Component fit 

Fig.QT]shows fits to the azimuthal velocity distribution (green data 
points) for the the four different subsamples and distance calibra- 
tions we use. Both columns show the same data, the left column 
on a linear scale, the right column on a logarithmic scale to exam- 
ine both centres and wings of the distributions. We apply a simple 
Gaussian halo and a single non-Gaussian disc component. All fits 
were done for -250 km s < V < 280 km s -1 , a limit that stays out 
of the regions dominated by noise, but reaches on the retrograde 



side still well into the region where any suspected retrograde halo 
component would be influential. To account for the observational 
errors, we folded the underlying distributions with a sum of Gaus- 
sians (with a spacing of lkms -1 ), of which the relative weights 
were derived globally for each component from error propagation 
on the single objects (stars with V > 50 km s _1 attributed to the disc, 
objects with V < lOkms -1 to the halo) concerning proper motion 
and radial velocity errors. This treatment is a bit crude and we sus- 
pect that it underestimates errors on the left wing of the distribution 
and overestimates them on the right wing, which we identify as the 
most likely reason for the overshooting of the model against the 
data at V > 270kms~ 1 . In a second step we folded with another 
Gaussian magnitude error of 0.25 mag to account for uncertainties 
in the intrinsic brightness and thus distance of the stars and there- 
after a third Gaussian of 7 km s -1 to account for velocity crossovers 
via the distance uncertainties. 

For the disc component we make use of the analytic formula 
of Schonrich & Binney (in prep.). The underlying assumption is 
that of an isothermal sheet through the Galaxy with increasing ve- 
locity dispersions towards the centre, i.e. the stellar populations at 
each galactocentric radius are given a specific vertical and horizon- 
tal energy dispersion. For those populations their likelihood to be 
in the solar annulus and their local scaleheight can be estimated as- 
suming a simple potential with constant rotation speed. We used a 
solar galactocentric radius Rq = 8.0 kpc, a circular rotation speed of 
v c = 220kms~', a disc scalelength of R d = 2.5 kpc, a scalelength 
for the vertical dispersion of R a „ = 5.0 kpc and for the horizontal 
dispersion as R a = 7.5 kpc. The adiabatic correction index (index 
dependent on the shape of the potential that describes the change of 
vertical energy along orbits that are extended over different galac- 
tocentric radii) was set to 0.5. As we are using the adiabatic cor- 
rection without a recorrection for energy conservation, the disc lo- 
cal scale height gets moderately underestimated, an effect that we 
partially cope with by setting the adiabatic correction index to 0.5, 
slightly below the expected value for the upper disc. To simplify the 
calculation we summed up disc contributions at equal parts at al- 
titudes z = 400, 700, 900, 1 100, 1400, 1 800, 2400 pc. The Gaussian 
halo component was set at rest (V = 0). The five free parameters of 
the fit were disc and halo normalization (Ni, = 530, N d = 16691 for 
the entire C10 sample, N d = 14403, N h = 312 for the C10 dwarfs, 
N d = 13031, N h = 262 for the adopted main sequence calibration, 
N d = 13254, A';, = 255 for the age-dependent calibration), halo 
azimuthal velocity dispersion (84.7, 80.9, 65.4, 68.9 km s~'), disc 
local horizontal dispersion (cr = 43.2,41.4,41.5,40.1 kms~'), 
and local scaleheight (ho = 571, 499, 466, 439 km s~'). As seen in 
Fig. [TT] apart from a weakness on the high velocity side, which 
presumably derives from the use of a global observational scat- 
ter that could overestimate the uncertainty for disc stars, the two- 
component approximation gives decent fits (red lines) to the veloc- 
ity distributions for the dwarf samples. The relative halo normal- 
ization varies between the different datasets as the most metal-poor 
stars are more likely to be classified into the turn-off or subgiant 
bands. 

Of most interest for this discussion is the shape of the (retro- 
grade) halo velocity distribution. For the entire C10 distance sam- 
ple (left column) we can clearly identify the bump that encouraged 
C07 and C10 to fit a separate velocity peak starting from around 
-200 km s" 1 . This anomaly becomes especially apparent in the log- 
arithmic plot and is confirmed by Poisson loglikelihood values rul- 
ing out equality of the theoretical (simple Gaussian halo) and ob- 
served distributions at highest significance. 

We point out that even this strong anomaly in the fit alone 
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Figure 10. Metallicity distributions at different values for kinetic energy in the different distance prescriptions. From top to bottom we show the entire CIO 
sample, the dwarf stars from CIO, the adopted main sequence calibration and the Ivezic et al.(2008) age-dependent calibration. The left column displays a 
separation by the entire kinetic energy, while the centre column makes only use of the vertical velocity W. The right column shows the distributions using the 
vertical energy part for retrograde stars (V < -10kms~') to reduce the disc contamination. Velocities are taken in kms~'. Error bars indicate the Poisson 
errors. 



would not be a sufficient justification for a second physical com- 
ponent, especially not for a retrograde component. There is no 
reason to firmly believe that the Milky Way halo or its possible 
components should have a strictly Gaussian velocity structure; the 
disc certainly can not be adequately described by G aussian fits 
(a discussion of this can be found in IStr6mberd[l927l) and simi- 
larly the halo might have a more complex velocity distribution. The 
only way to cleanly identify a retrograde velocity distortion on the 
Galactic halo from kinematics would be to prove a difference to the 



prograde halo tail. As this is, however, impossible at the required 
accuracy due to the disc contamination, we conclude that without a 
clean halo-disc separation any attempt to use the azimuthal velocity 
distribution for proving a retrograde halo component is unreliable. 

Besides its doubtful significance for indicating a separate 
component, the extended retrograde tail anyway diminishes when 
we turn to the dwarf samples, which do not contain the at least 
partly unphysical turn-off stars: restricting the CIO sample to their 
dwarf stars (centre column of Fig.ll U the bump below -200 km s~ 1 
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Figure 11. Velocity distributions for stars (light blue, with Poisson errors) in the calibration sample for different selections and distance prescriptions: From top 
to bottom: The entire C10 sample, the C10 dwarf stars with log(g) > 4.0, the dwarfs with log(g) > 4.1 using the adopted main sequence calibration distances 
and the age-dependent formula from Ivezifi Z et al. (2008). To show both distribution centre and wings we contrast the linear scale (left) with logarithmic plots 
of the distributions (right). All distributions are fitted by a simple model (red line) with a Gaussian halo and a single (non-Gaussian) disc component. We did 
an error propagation from the values given in DR7 which is - together with a magnitude error of 0.25 mag - folded onto the original distribution (light blue). 
The error broadens the peaks and thus lowers the maximum count rates especially of the disc peak. The dark blue line shows the underlying Gaussian halo 
component and the purple curve the analytic disc component. 
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disappears. There is a weak surplus of stars between -220kms~' 
and -280km s - ', which could, however, just be a statistical fluctu- 
ation. This is confirmed by statistics, which have the p-level (mea- 
sured under the assumption of Poisson errors) for equality of the 
fit and the measured distribution at 0.55 for -250kms~' < V < 
20km s _I (and 0.20 when we extend the range to -300kms~' < 

V < 20kms~') meaning that when drawing realisations from the 
given theoretical distribution, about 55% of the samples would 
be more different from the theoretical distribution than the cur- 
rent one. Regarding our imperfect treatment of distance errors the 
quality of the fit is rather surprising. Using the main sequence dis- 
tance calibration the distribution gets even more contracted due to 
the shorter distances. Apart from the slight surplus of stars left of 

V ~ -320 km s -1 , which just appear to be some sample contam- 
ination (indeed a distance test on the 100 objects with lowest V 
velocities still reveals a 2a distance overestimate) and the single 
bin with 8 stars around V = -260kms~' no anomaly is trace- 
able any more, confirmed by statis tics. Using the age-dependent 
calibration from llvezic et al] d2008T ) the picture is similar. The fit 

looks even better down to V 300 km s~' beyond which there is 

a weak surplus of about a dozen objects between V = -400km s -1 
and V = -300kms~'. This can be attributed to misassignement 
of stars between the catalogues for proper motions and a probable 
non-Gaussian far tail of the proper motion error distribution. 

It can also be argued that an attempt to fit the disc velocity 
distribution via a Gaussian fit increase the need for a spurious sec- 
ond halo component, especially as the large (unavoidable) veloc- 
ity errors limit the information on the real shapes of the underly- 
ing velocity distributions: as the disc velocity distribution is nat- 
urally skewed towards lower velocities, a Gaussian fit will drop 
too steeply towards low rotation velocities. This invokes a second 
Gaussian for the disc (thin-thick discs), which will in most cases 
show the same problem again, forcing the halo component a little 
bit up into the prograde regime. This again increases the need for 
the creation of an artificial retrograde halo component. In this case 
we avoided this problem by using a more physical fitting formula 
for the Galactic disc, yet there was still noticeable influence by the 
disc fit onto the parameters of the simple Gaussian halo component. 

Overall we can state that any striking excess of highly ret- 
rograde stars disappears from the distribution when we limit the 
sample to the more reliable dwarf stars. There is a good agreement 
between fits and data for the dwarf samples, apart from a misfit at 
high velocities and a slight surplus at the lowest velocities against 
a Gaussian fit, which, however, rather looks like imperfect error 
handling or a probable contamination of the sample with misiden- 
tifications in between the different catalogues delivering the proper 
motions or a non-Gaussianity in their contamination. Replacing the 
questionable Gaussian analysis for the disc by a more physically 
motivated formula and by applying an error propagation we have 
shown that neither a second halo nor a second or even third disc 
component are required for explaining the azimuthal velocity dis- 
tribution. 

5.3 Velocity dispersions 

Using the adopted main sequence approximation and restricting the 
sample to 2020 retrograde stars (V < -10kms~') we can give a 
tentative estimate for the halo kinematics of cry ~ 157 ± 10 km s -1 
and o"iy ~ 75 ± 8 km s~ 1 . The latter is close to the azimuthal velocity 
dispersion of a v ~ 70kms~' from the last section. In these values 
we accounted for the errors reported in the SEGUE pipeline plus 
a dispersion of roughly estimated 15% in the distances. These act 



to reduce the derived velocity dispersions from the actually mea- 
sured ones. Those values are considerably smaller in V and W than 
what was given by C10 partly due to the spurious crossover into 
W velocities having diminis hed. We point out t hat these values 
agree well with the results o f Kep lev et alj (2007) and lSmith et al.1 
(2009). Although w e use the same sample with the same distance 
scale as ISmithetal.1 120091) we cannot confirm their quoted errors of 
2 k ms" 1 , since the sys tematic uncertainties by the distances within 
the llvezic et"aT] (j2008) method and sample cleaning are currently 
too large to assess dispersions on a scale better than 5%. 



6 CONCLUSIONS 

We have described how errors in distance estimates result in an ap- 
parent systematic retrograde motion of the Galactic halo, an effect 
to which the SEGUE/SDSS sample is especially prone by its strong 
poleward orientation. The general problem of distance biases simi- 
larly applies to any study that makes use of pro per motion-based es- 
timates. We find that the distance derivation of Carollo et al. 
and lCarollo et al.l d2010h is flawed by sorting stars into unphysical 
positions in the HR diagram: objects are placed between the sub- 
giant and dwarf sequences in positions that would require stellar 
ages in excess of the age of the universe. Despite the elegance of the 
general idea to sort stars into known sequences according to their 
estimated gravities, the method itself and the used gravity cuts are 
not well supported by measurements. Moreover there is no "turn- 
ofF'-sequence, but turn-off stars are populating a region that spans 
of order 1 magnitude in luminosity. In this light the statement by 
C10 to have distances precise to about 10 - 20% is an unsupported 
claim. 

From the distances kindly provided by Carollo et al. we calcu- 
lated back to their assumed absolute magnitudes and found sys- 
tematic differences of ~ 0.2 to 0.3 mag and a large scatter for 
metal-poor main sequence stars towards the adopted main sequence 
calibration as wel l as towards the age-dependent calibration by 
Ivezi c et al. I < l2008l) . also far to the red side of the suspected turn- 
off region. The adopted main sequence calibration is only slightly 
fainter than the theoretical BASTI isochrones. 

We have shown in Section POl that the claim by C07 and C10 
to have found a counter-rotating extended tail of the halo is largely 
caused by unphysical assumptions about locations of stars in the 
HR diagram, by magnitude uncertainties in the turn-off stars and 
by the use of a too bright main sequence calibration. The correct- 
ness of this tail can be ruled out by statistical tests, as described 
in Section|4] The tail diminishes when we limit the C10 sample to 
dwarf stars and disappears when we make use of the better founded 
llvezic et"aT1 d2008l) calibrations, which are consistent with fiducials 
and isochrones. In Section l5?2l we demonstrated that the halo dis- 
tribution for the dwarf samples regardless of the applied distance 
determination can be fit by a simple Gaussian component. 

We have also shown that in the DR7 pipeline stars with lower 
metallicities are shifted towards lower gravities, considerably in- 
creasing their fraction among the thought-to-be turn-off stars. Fur- 
the r the magnitude diff erence for their main sequence stars against 
the llvezige t al. (2008) main sequence calibration grows towards 
lower metallicities. The stronger prevalence of distance errors at 
the metal-poor end of the metallicity distribution will thus give any 
spurious counter-rotating tail members a biased metallicity distri- 
bution. 

Finally we have shown that the claim of C10 that the counter- 
rotating component members reach to higher altitudes can as well 
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be traced back to distance determinations: the dispersion of the ver- 
tical velocity component is significantly increased by their distance 
errors, though due to the polewards sample orientation this effect 
would at first order be smaller than for the other velocities. As 
shown in Section 10 simultaneously the W velocities of the halo 
stars with distance overestimates are artificially increased by of or- 
der 50 km s via spurious velocity cross-over terms from the helio- 
centric azimuthal velocities in the derivation. The effect is strongest 
for the most strongly retrograde objects (they have the largest he- 
liocentric Vh) and is aggravated by their selection for stars in galac- 
tocentric radius 7 kpc < R < 10 kpc. This colludes with their metal - 
licity dependent distance bias (see above) to produce their findings 
of decreasing metallicities at high altitudes. There is a slight excess 
of more metal-poor stars in a single velocity bin at high vertical ve- 
locities, which is not mirrored by the behaviour at high total kinetic 
energies. We argue that this is most likely a r eflection of a well - 
kno wn local stream tha t has been identified bv lHelmi et al.l J 19991) 
and lKeplev etai] J2007h . 

Another source of error is the modelling of especially the 
Galactic disc azimutha l velocity d i strib ution by Gaussian compo- 
nents. It was shown by Stromberg] J 19271) that Gaussian modelling 
of the Galactic disc lead to unphysical results and the identification 
of spurious components on the low rotation side because of the ex- 
tented tail. As the skewed V velocity distribution enforces in most 
cases the introduction of a second Gaussian component, that - being 
a mere artifact by wrong assumptions - can then be misinterpreted 
as physical reality, Gaussian modelling of the Galactic disc in a 
combined disc and halo sample can wrongly force the halo com- 
ponent into the prograde regime to compensate for the two steeply 
falling disc terms. Consequently this then creates the need for in- 
ference of a retrograde component to compensate for the bias. 

We also argue that magnitude based distance assessment 
schemes introduce a velocity bias that resembles the Lutz-Kelker 
bias: If the error in absolute magnitudes follows a Gaussian distri- 
bution, the distance error distribution will thus form an extended 
tail that grows stronger with increasing dispersions. Via the proper 
motion part in the determination of space velocities, which is pro- 
portional to the estimated distance, measured velocities develop ex- 
tended tails away from the solar motion. For the V velocity dis- 
tribution of especially the halo this gives the halo an asymmetric 
velocity distribution with a longer tail in the retrograde regime, a 
process that can expl ain the moderate asymmetries found e.g. by 
iNorris & RvanUl989l) . 

Finally we note that it is by no means imperative that the halo 
have a Gaussian velocity distribution. In this light it is rather sur- 
prising that our simple Gaussian halo component can fit the data 
so well. Even if there were deviations from Gaussianity this would 
alone be no convincing sign for a separate component. A more con- 
vincing indication would be a proven difference between the pro- 
grade and the retrograde tail of the halo azimuthal velocity distribu- 
tion, but disentangling this from disc contamination on the prograde 
side will be difficult. 

Recently Carollo and collaborators submitted a rebuttal pa- 
per jBeers et al]|201 ll) claiming that our analysis presented here be 
wrong. Instead of discussing all our arguments their revised anal- 
ysis relies on two central claims: They state that the distance scale 
adopted by us is wrong and that there is an asymmetric halo az- 
imuthal velocity distribution for their metal-poor stars, neither of 
which we concur with. 

Concerning the distance issue we stress that our conclusions 
are va lid for the Carollo et al. (2007, 2010) and both Ivezic et al. 
j2008h distance calibrations; our work does not rely on a single 



distance scale in contrast to claims made by iBeers et al] d201 lb . 
Beers et al. c riticize us for adop ting the incorrect main sequence 
calibration o f llvezicetal.ld2008l) but failed to note that we actually 
stretched this calibration in the same direction of their preferred 
one by increasing the luminosities by 0. 1 magnitudes and account- 
ing for alpha-enhancement by increasing the measured metallicities 
by 0.2 d ex (Sect. 2). I mport antly, we have also made use of their 
preferred llvezic et alj d2008l) calibration (here denoted IvzA7) and 
find no significant differences (e.g. Fig. 1 1). Finally, we have made 
use of direct isochrone distances, which fully corroborates our find- 
ings. 

As for the azimuthal velocity distribution of the most metal- 
poor stars, we re-emphasize that any magnitude-based distance 
scheme invokes a bias in the inferred distances and thus an asym- 
metric azimuthal velocity distribution by definition; this effect is 
akin to the well-known Lutz-Kelker (1973) bias and is illustrated 
in Fig.QT] In view of a large magnitude scatter (which their metal- 
poor stars clearly have, see their Fig. 5) a Gaussian fit is inappropri- 
ate due to the missing error handling. In this light it is not surprising 
that their new revised parameters are quite different from their orig- 
inal results (e.g. for z„, av > 5 kpc their outer halo mean velocity rose 
from -128 km s _I (cf. Table 1 in C10) to -59 ± 20 km/s). Finally, 
we argue that moving a considerable fraction of the wrongly iden- 
tified turnoff stars up to the subgiant/giant branch as done by Beers 
et al. will make the distance overestimate for misidentified dwarfs 
among them even more severe. 

In summary our criticism of the C07, C10 works remains in 
full: our in-depth re-analysis of their data with different distance 
calibrations and a proper error handling reveal no convincing evi- 
dence for a dual halo. 

The systematic distance uncertainties make it dangerous to 
draw a definitive conclusion for the strength or existence of a pos- 
sible counter-rotating halo component. All current distance calibra- 
tions have problems and need improvement before a method along 
the lines used by C07 and C10 can be attempted. We would like to 
stress that we do not and would not want to rely on either of them. 
Two central conclusions can, however, be drawn without having 
to trust any of the different distance calibrations: Even on the C10 
or C07 sample using their distances, no reliable detection of any 
non-Gaussianity in the halo, be it a counter-rotating halo or not, is 
possible on the examined data set in any of the dwarf samples. If a 
separate component gets detected on a larger sample in the future, 
it should be significantly weaker than what was claimed by C10. 
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8 APPENDIX 

8.1 Kinematics and geometry 

We use up-to-date values for the basic constants of our Galaxy. 
The rotation speed is assume d to be 220 km s -1 in concor- 
dance with recent results of iKoposov et alj {2010). We apply 
the recent determinati ons of solar mo t ion re lative to the local 
standard of rest from ISchonrich etafl bold) , which is v n = 
(H.1^75, 12.24^, 7.25!jJ;^)kms- 1 with additional systematic 
uncertainties of ~ (1, 2, 0.5) kms~'. Neither the Galactic ro- 
tation rate nor the solar galactocentric radius R Go are very 
well determined, but the angular motion of Sagittarius A* is 
jReid & BrunthaleJ |2004|) . The galactocentric radius of Rq.o — 
8.5 kpc assumed in C10 is inconsistent with their assumption for 
the rotation speed. We adopt Rg.q = 8.0 kpc, which is in concor- 
dance with most measurements, and coincides with the most recent 
trigonometric parallax determination for Sagittarius B2 jReid et all 
120091) . 



8.2 Distance calibrations and metallicities 

For adopted dis t ances in SDSS colours two alternatives exist: 
The llvezic et alj (2008) calibration (see Appendix therein) or the 
Beers (2000) calibration. A third possibility would be using the 
isochrones directly, which would allow for a statistical implemen- 
tation of the subgiants and also allow for natural shape changes: 
however, their ha ndling is beyond the scope of this work. The 
llvezic et "all (2008) main sequence calibration uses the formula: 

M r (g - r, [Fe/H]) = (3) 

= 1.65 + 6.29(g - Oo - 2.30(£ - i% - l.ll[Fe/H] - 0.1 8 [Fe/H] 2 (4) 

where M r is the adopted absolute magnitude and (g - Oo is the 
dereddened SDSS colour index. A short assessment of this formula 
reveals that it fits the zero age main sequence relatively well in 
the required colour range at low metallicity. Apart from the short- 
come that the subgiants are not considered, the metallicity depen- 
dence is not well matched at the high metallicity end. Further the 
colour dependence is not perfect at the level of precision required 
for distances and there are no crossterms between metallicity and 
colour, i.e. there is no implementation of metallicity-dependent 
shape changes. Thus this relation has to be used wit h caution. A 
compa rison with the BASTI isochrones reveals that the Ivezic et al. 
( 120081) calibration at low metallicities is slightly fainter than the 
isochrones without alpha enhancement. To correct for alpha en- 
hancement we slightly increase the metallicity of all metal-poor 
stars (with [Fe/H] DR1 < -1.0) by 0.2 dex and let this correction go 
linearly to zero between [Fe/H] DR j = -1.0 and [Fe/H] DR ~i = 0. 
In addition we increase all luminosities by 0. 1 mag to reconcile 
the calibration better with the main s equence of M92 accord- 
ing to the fiducials by lAn et alj d2008l) . to avoid any suspicion 
to have buried the counter-rotating tail by a too faint calibration, 
and as we detected a minor distance underestimate in our statis- 
tics without this shift. Further the calibration is now just slightly 
brighter than the isochrones as seen in Fig. [4] which is a desired 
effect at those metallicities, as there are indications that metal-poor 
isochrones could underestimate the luminosities (or equivalently 
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Figure 12. Vertical dispersions vs. metallicity revisited. The lines are the 
same as in Fig. [9] Additionally we plot in black the same derived via direct 
isochrone determination of distances. 



overe stimate the effective te mperatures) of lower main sequence 
stars dCasagrande et alj|2007l) . 

To account for age effects frvezic et "all d2008h also suggested 
an "age-dependent" calibration, which rises more steeply towards 
the blue side (Eq. A 7 in their appendix). To meet concerns of our 
referee we show all relevant statistics also in the light of this other 
calibration. This calibration is steeper in colour than the main se- 
quence calibration, i.e. it is brighter at blue colours and fainter at 
red colours, following an intention to cope with the steepening of 
isochrones near the turn-off. However, the steepening of isochrones 
only applies near the turn-off, while this relation is globally inclined 
against the main sequence of the isochrones as well as against the 
main sequence calibration. As this relation thus does not follow the 
blueward shift of the turn-off towards lower metallicity, a relative 
overestimate is expected for the distances of the blue stars at low- 
est metallicity. Vice versa the distances to metal rich turn-off stars 
might be underestimated. 

We tested all our results for the CIO dwarf sample and also 
for cutting in colour and different gravity selections to delineate 
the impact by the turn-off region. As can be seen from the central 
panel in Fig. [4] the difference between the adopted main sequence 
calibration and the CIO magnitudes persists al so to the red side of 
the turn-off region and is also present for the llvezic et al.l ( 2008) 
age-dependent relation. 

As an even fourt h distance determination we made use of 
the BASTI isochrones jPietrinfemi et alj|2004 l2006h . To accom- 
plish that we account again for alpha enhancement with the same 
prescriptions as for our adopted main sequence relation. For each 
star we choose the closest 12.5 Gyr isochrone in metallicity from a 
d ense grid kindly provided by S. Cassisi for our age determinations 
in lCasagrande et al]d201 lh . On this isochrone we choose the r-band 
magnitude from (g-i) on the main sequence. When the turn-off is to 
the red of the stellar position, we extrapol ate from the turn-off point 
of the isochrone following the shape of the llvezic et alj ( l2008h main 
sequence relation. So these stars are placed at an extrapolated turn- 
off point. As one can easily see from Fig.[l2]in the central statistics 
there are no significant changes towards the other main sequence 
calibrations. We confirmed that no real difference was detectable. 



CIO used a slightly different metallicity scale from DR7, 
which was reasoned to balance out a possible metallicity overes- 
timate on the lowest metallicity end: 

[Fe/H]car = "0.186 + 0.765 [Fe/H] DR7 - 0.068[Fe/H] DR1 2 (5) 

We argue that their claimed overestimate at low metallicities most 
likely reflects the impossibility of the pipeline to cope with the faint 
metal lines in this region, getting lost in the low resolution and low 
signal-to-noise ratios. This effect is quite similar to the loss of ac- 
curacy e.g. in Stromgren photometric metallicities. The formula 
has the dissatisfactory effect of aggravating the well-known bias 
of the DR7 pipeline to underestimate the metal content of metal- 
rich stars. This can be seen by comparison of local stars from the 
Geneva-Copenhagen Survey (Nordstr om et al. 120041) with high ver- 
tical energy, to the metallicity distribution of the SEGUE disc stars. 
Similarly the ugriz metallicity calibration fails to reproduce metal- 
licities al ready slightly below sol ar metallicity as can be seen from 
Fig. 1 1 injArnadottir et al. ( 2010). The latter problem is not of ma- 
jor importance for the halo, but applying the formula by Carollo et 
al. would exacerbate metallicity-induced errors on the disc popu- 
lation. We estimate that distances are only weakly affected by this 
correction on the low metallicity side (i.e. mostly for halo stars) as 
the sensitivity of stellar atmospheres t o the logarithmic metalli city 
scale gets lower. This can be seen in 



VandenBerg et al.l d2010h or 



for comparison in fig. 13 of lCasagrande et a" 1 20ld) . 



We use the DR7 metllicities throughout the paper and do not 
apply the correction from l[5}, but checked that our findings are not 
significantly affected by switching the distance scale. 

8.3 Deriving space velocities 

The dataset contains information on magnitudes, colours, the dis- 
tances, stellar parameters, radial velocities and proper motions. We 
can thus derive the velocities in the solar/local coordinate system 
by: 



Uh = -dsin(l)l - dcos(/) s'm(b)b + cos(fe) cos(/)v|| 
V/, = d cos(/)/ - d sin(7) sin(b)b + cos(b) sin(Z)v|| 
W), = d cos(b)b + sin(fc)i'n 



(6) 



where /, b are Galactic longitude and latitude, /, b are the proper 
motion components in / and b and d is the assumed distance to the 
stars. Of special interest is the degree of support of the measure- 
ment by radial motions, which are independent of distance biases 
and have smaller errors. Neglecting the geometrical extension of 
the sample, the quantity of interest is thus the term connecting the 
azimuthal velocity in the heliocentric frame V/, to the line-of-sight 
velocity V|| : 



Tj V = |cos(fc) sin(7)| 



(7) 



The term rj v is 1 where the V;, velocity is measured directly from the 
line-of-sight velocity vy i.e. at b = and I = 90°, 270°, its contours 
on the sample are depicted in Fig.Q] 

Due to the extension of the sample to a radius of more than 
3 kpc projected on the Galactic plane a small angle approximation 
cannot be reliably taken. Using the Galactic rest frame velocities in 
heliocentric coordinate U a = Ui, + U B , V a = V;, + Vg + V , W a = 
W* + Hq (where V G is the Galactic rotation speed, V B , U and W Q 
are the components of solar motion relative to the local standard of 
rest) the correction from the local coordinate system can be done 
via the Galactic angle 



dsm(l) cos(b) 

a = arctan( ) 

Rqo - acos(/)cos(o) 



(8) 
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between the line sun-centre to the line star-centre: 

V = U a cos a - V a sin a 
V = V a cos a + U a sin a - V G (9) 

W= W a 

Throughout the paper blank letters U, V, W denote the corrected ve- 
locities in a Galactic reference frame, which are useful for assessing 
kinematics, I//,, V/,, Wi, denote the velocities in the rest frame of the 
Sun and heliocentric coordinates, which are the native setting for 
exploring kinematic biases. 

8.4 Sample cleaning 

The DR7 calibration sample contains 42841 spectra. The steps for 
cleaning the sample and the subsequent reduction of numbers are 
listed in Table [T] subsamples used in this work and the plots are 
set in bold font. However, a considerable fraction of these entries is 
double, i.e. stars whose spectra have been taken several times, some 
for observational reasons as to improve the signal to noise, some as 
they were both listed among the photometric calibration stars and 
the reddening calibration stars. Therefore it is necessary to clean 
DR7 samples by identifying any measurements that are within 1.5 
arcsec of each other in both right ascension and declination or are 
within an angular distance of 3 arcsec and have a g band magni- 
tude difference below 0. 1 mag. We chose the entry that gives full 
information on proper motions and the latest measurement in case 
both contain it. This leaves 33023 unique calibration stars in our 
DR7 sample. The CIO sample was already cleaned by them and 
has the 4kpc distance cut (according to their distances) applied, 
but not the cut in galactocentric radius (which we do not apply ei- 
ther). Hence the number of unique stars in the CIO sample is only 
slightly reduced when we demand a cross-match by stellar position 
on our DR7 table for being less then 2.5 arcsec apart both in right 
ascension and declination. There are some stars dropping out be- 
cause of a missing cross-match and there were 52 candidate double 
entries with identical position in the CIO sample. Among those ob- 
jects some are ~ 0.3 mag fainter than their second entry and the 
corresponding entry in DR7. In total we found 41 stars which are 
~ 0.3 mag fainter in apparent magnitude than their counterparts in 
DR7. We checked that none of these differences has a significant 
impact on the results. 

We also removed from the sample the stars with signal to noise 
ratio S/N < 10 and those which are flagged by the SSPP for spec- 
tral abnormalities, for colour mismatches or for being a suspected 
or proven white dwarf, as well as those with particularly unreliable 
radial velocity measurements or SSPP parameters. We thus cleaned 
the C10 sample from ~ 1700 flagged stars. Superficial tests did not 
show any obvious problems caused by those stars. Following C10, 
we also excluded all stars witho ut determined pro per motions and 
"clean" the sample according to Mu nn et al] d2004l) requiring that 
both cr RA and cr DEC < 350 mas, with an additional requirement of 
those stars to have quoted errors in each proper motion component 
of < 5masyr~'. We checked on the C10 sample that 224 objects 
in their sample not passing this cut did not cause any noticeable 
biases. The full cleaning gives a final DR7 sample of 28844 stars, 
which we will use when plotting "all" stars in DR7. In plots of 
the full C10 sample we use its counterpart of 21600 stars that al- 
ready fulfills at their distances the condition of having all distances 
smaller or equal to 4kpc. 



original sample entries 


42841 


— 


cleaned sample by C10 


— 


23647 


unique & cross-matched 


33023 


23553 


unflagged 


29655 


21828 


S/N > 10 


29638 


21825 


4500 < T cff < 7000 


29601 


21823 


[Fe/H] DR7 ,v n fine 


29584 


21819 


prop, motion fine 


28844° 


21600 


log(g) > 4.0 


17365 


15023 


log(g)>4.1 


13880 


12120 


d h , z < 4kpc,log(g) > 4.0 


15808 


14763 


d lvl < 4kpc,log(g) >4.1 


12678 6 


11894 


d Car < 4 kpc 




2160IF 


d Car <4kpc,log(g)>4.0 




15023'' 



Table 1. Numbers of stars at different cuts in the two samples. The quality 
cuts are applied cumulatively from the first row until the horizontal line. 
Below the horizontal line we show selected subsamples with gravity and 
distance cuts. The cuts down to the quality cut in proper motion are applied 
successively. Below the horizontal line we show the effects of differnt dis- 
tance estimations and cuts on the number of remaining stars. 
Subsamples used in the paper: " "all star sample", * "Ivezic dwarfs", c Car- 
ollo all, d Carollo dwarfs. 



